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Simultaneous optimization of multiple quantum objectives is often considered a demanding task. 
However, a special circumstance arises when a primary objective is pitted against a set of secondary 
objectives, which we show leads to invariant behavior of the secondary objectives upon the primary 
one approaching its optimal value. Still, practical relationships among the objectives will gener¬ 
ally lead to a threshold, beyond which system re-engineering is required to further increase the 
primary objective. This finding is of broad significance for reaching high performance in quantum 
technologies. 


The control of complex quantum phenomena [1] 
can entail simultaneous optimization of multiple ob¬ 
jectives in many applications. Nevertheless, a com¬ 
mon circumstance is to seek maximization of some 
primary measure J to high fidelity. For example, 
in the emerging field of quantum technologies [2] 
J may be defined as the success of a quantum in¬ 
formation processing gate [3], or a measurement 
from a metrology or a quantum sensing experiment 
[4, 5] . Additionally, while seeking high fidelity for J, 
there will typically be additional secondary objec¬ 
tives {AT,} such as performing control with the least 
amount of energy, preserving coherence, or avoid¬ 
ing unwanted noise effects. Optimization of the pri¬ 
mary and secondary objectives are often performed 
by seeking tailored, applied fields, and consideration 
of J and {Ki} objectives will generally involve a 
tradeoff. However, this letter will point out an at¬ 
tractive feature of unitary quantum control impart¬ 
ing a special general character to this tradeoff in the 
desired high fidelity regime for J. 

We assume that the physical circumstances de¬ 
mand that J be nearly optimal, and we aim to 
discern the tradeoffs when additionally optimizing 
{Ki}. Characterizing such tradeoffs is commonly 
called the Pareto problem, where typically increase 
of one objective value through changing the controls 
cannot be performed without penalty to the other 
objectives [6, 7]. Such circumstances pose a difficult 
situation, as no operating condition may be satisfac¬ 
tory for meeting all concerns [8-12]. 

Here we demonstrate that rather than finding a 
tradeoff, at high values of J an invariance settles in 
for the value of all secondary objectives {Ki}. This 
behavior holds for any well posed primary quantum 
objective J, and is important in the growing efforts 
to develop quantum technologies. To illustrate, we 
will demonstrate its implications in three physical 


circumstances (see eqs. 4-6) where high fidelity of 
some operation J is necessary. Spin systems will be 
used in this regard. As the two objectives we take 
the fidelity of J and a measure of robustness Kp to 
noise disturbances of type /? in the system. Numer¬ 
ical calculations will show that for increasingly high 
fidelity of J, the robustness Kp to either detuning 
noise or control noise freezes in at an invariant value. 
The same conclusions would apply to analogous op¬ 
timization operation directly performed in the labo¬ 
ratory. In general, the invariance of Kp to increasing 
fidelity of J will reach a threshold, whereby proceed¬ 
ing further with J will call for improving Kp that can 
only be achieved by appropriate re-engineering of the 
system. Depending on the circumstances, dealing 
with infringing {Ki} could call for re-engineering of 
the control sources, the physical system, the detec¬ 
tors, etc. While this invariant behavior is indepen¬ 
dent of search technique, it does provide algorithmic 
advantages, as even a modest ensemble of Monte- 
Carlo (MC) samples of nearly optimal controls can 
readily identify the best invariant value of I\p. 

The invariant behavior is shown to arise due to 
the generally expected, unencumbered convergence 
toward the optimal value of J along a path on the 
quantum control fidelity landscape , defined as the 
functional mapping from the controls e(t) to the ob¬ 
jective J. As background, the success of a growing 
number of optimal control experiments led to the 
formulation of a key theorem, referred to here as the 
“landscape principle,” demonstrating that the topol¬ 
ogy of well-conditioned quantum control landscapes 
allows for facile determination of optimal controls 
[13-15]. This principle provides a foundation to ex¬ 
plain the growing successes of quantum control [16], 
and at this juncture very few exceptions are known 
to this principle, thereby providing a solid basis for 
its practical exploitation [17]. 


The convergence of an optimal control e(t) to pro¬ 
duce a high value of J is observed by parametriz¬ 
ing the control by a search variable s > 0. In this 
fashion, traversal over the landscape will start from 
s = 0, where J is generally at low fidelity, out to a 
final value Sf = oo (i.e., in practice some finite value 
of Sf where J is near optimal, taken as J = 1). A 
gradient-based optimization along s will always seek 
to increase J through 


dJ__ f T 6J de(s,t 
ds J Q Se(s,t ) ds 

which is guaranteed by [18, 19] 

de(s,t) 5J 
ds Se(s,t) 


(1) 

( 2 ) 


Numerous simulations and experiments have demon¬ 
strated that the landscape principle enforces that 
norm of the gradient SJ/5e(s,t) has a Gaussian-like 
shape as a function of s that eventually decreases as 
s — y Sf until a critical point at the landscape top is 
approached where 5J/5e(sf,t) ~ 0, Vt. The condi¬ 
tion in eq. (2) implies that near optimal fidelity, 


lim 

S—tSf 


de(s, t) 
ds 


— 0 , 


( 3 ) 


where ever decreasing variations in e(s,t) occur 
for s —> Sf. Consequently, ancillary objectives 
Ki[e(s,t)] that solely depend upon the controls also 
converge at high s values. 

While traditional Pareto “front” tradeoffs will 
likely exist when initially climbing the control land¬ 
scape, the asymptotic behavior of all individual 
search trajectories via eq. (3) (i.e., each guaranteed 
to reach high fidelity for J by the landscape prin¬ 
ciple) likewise assures a high degree of invariance 
between the fidelity of J and secondary objectives. 
Figure 1 demonstrates this global behavior by show¬ 
ing how a MC sampling of control trajectories climb¬ 
ing the landscape can readily reveal convergence to 
the highest invariant value of field noise robustness 
K e at ever lower fidelity error Ej (i.e., higher J) 
given as Ej = 1 — J. Maximization of the secondary 
objective K e is seen as a single, flat curve labeled + 
in the regime of high fidelity (i.e., log (Ej) < —6) 
where K e plateaus [20]. 

However, this invariance hides a competition, as 
a threshold value of J will occur where Ej = — Kp 
(Fig. 1, black dots for the highest value of Kp) , 
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FIG. 1: Color online. Field noise robustness 
measure in eq. (8) for representative collection of 
Monte Carlo optimization trajectories. Each curve 
arises from a single run up the Pi _>.2 landscape 
with a randomly chosen initial field. Convergence 
to high J (i.e. at smaller error Ej) demonstrates 
an invariance, evident in each line culminating in a 
single, non-dominated and increasingly flat 
trajectory for K t (labeled with + signs). Threshold 
point where Ej = —K e is given as a dot. The 
improvement of the threshold point is shown here 
by re-engineering noise statistics to increase K e 
(labeled with o signs). Inset, the probability 
distribution p (K t ) of 2000 MC samples at 
log(Ej) = —7.5, showing a preference to sample 
the largest invariant K e value. 


such that pushing further on fidelity can only be 
achieved by some suitable re-engineering of the sys¬ 
tem to increase Kp. Such an effect is demonstrated 
by the improved trajectory in Fig. 1 (o signs), in 
which the noise statistics were re-engineered (i.e., a 
longer noise correlation time), leading to an increase 
in Kp and an accompanying lowering of the thresh¬ 
old. Similar threshold behavior on any secondary 
objective likely occurs in all quantum control prob¬ 
lems (see Supplementary Material). 

We consider three types of primary objectives: 

p i^f = l(/l^( T )l*)| 2 State ( 4 ) 

(O) = Tr[p(T)0] observable (5) 

F w = 1-^-\\U(T)-W\\ 2 unitary (6) 

J is taken as the objective in either eqs. (4), (5) 
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or (6), which range here from worst fidelity (J = 0) 
to optimal (all normalized to J = 1). Control of 
any iV-level closed quantum system manipulates the 
dynamics of a unitary propagator U(t) = U{t,0), 
and correspondingly either its state \ip{t)) = U(t)\i) 
or density matrix p(t) = U{t)piW(t). U(t) is 
driven by a Hamiltonian H(t) through solution of 
the Schrodinger equation with short time steps At, 

U(t + At) = (7) 

For our numerical demonstrations, M coupled spins 
form the Hamiltonian. Each spin has a field-free 
drift term with transition energy w* between the spin 
states |0) = [1,0] T and |1) = [0,1] T , a control field 
£i(t), and collectively isotropic Heisenberg coupling 
terms of strength Ji j linking the spins (dynamical 
details provided in Supplementary Material). 

Maximization of the fidelity for two and four-level 
systems ( i.e. M = 1 or 2, respectively) are consid¬ 
ered for the following objectives: Pi-> 2 , 4 , (cr x ) 

{Pi = |0)(0|), and (ai 1) ) ( p % = |00)(00| ), where (a x ) 
and (cri 1 '*) are all rescaled to range from 0 to 1, with 
1 being optimal. For unitary transformation control 
in eq. (6) , we consider one-qubit Hadamard (Fh) 
and two-qubit CNOT gates (Fcnot)- Two level 
dynamics are simulated using eq. (7) over the in¬ 
terval T = 1 (At = 0.01), while four level dynamics 
are performed over T = 6 (At = 0.06). wi = 20, 
W 2 = 24, and a weak spin coupling is chosen as 
Ji 2 = 0.2. Optimization is performed via the D- 
MORPH gradient algorithm by integrating eq. (2) 
with a fourth-order Runge Kutta integrator (Mat- 
lab’s ode45) over a suitably long s interval [18, 19]. 

A measure for robustness to noise that is valid 
at high fidelity for the primary objective J was re¬ 
cently examined in [21], defined as a convolution of 
the Hessian 'Hp{t,t') = 5 2 J/5f3{t)5(3{t') of J and a 
noise correlation function R = {5/3(t)6(3{t')) p char¬ 
acterizing a stochastic parameter in the Hamilto¬ 
nian f}(t) (/3 is either e, or w,) subject to noise as 
/? (t) —» /3{t) + 5/3{t). The robustness measure Kp is 
the expected change in J due to noise in /3, 

K P = ^ W-p,R) 

= i f f 'Hp(t,t')R(t,t , )dtdt'. (8) 
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Kp values here are negative at high J fidelity, and 
robustness is enhanced when Kp is maximized to¬ 
wards zero value. 


While the degree of robustness depends upon the 
form of R [21], which is often not well known a pri¬ 
ori , the link to the landscape principle guarantees 
that the near invariance behavior of Kp should not 
significantly depend on the choice of R. For demon¬ 
stration purposes we restrict ourselves to Wide Sense 
Stationary (WSS) processes whose correlation func¬ 
tion decays in time delay t = t — t', as these well 
represent common noise forms found in quantum 
systems, such as 1//, l// 2 and white noise [22- 
25]. A decaying exponential noise correlation func¬ 
tion R(t,t’) = A 2 exp (— \t — t'\/a) is used, with a 
correlation time of a = 1 utilized to represent mid¬ 
frequency noise (a = 2 in Fig. 1, o curve) [21], and 
the noise strength is arbitrarily chosen as A 2 = 1(R 4 . 
For the four-level system, we assume similar but un¬ 
correlated noise statistics from each spin such that 
robustness is quantified as a sum of contributions 
from each field or energy level as Kp = Kp 1 + Kp 2 . 

Calculations were performed that seek to mini¬ 
mize the error Ej in the primary objective and max¬ 
imize the robustness measure Kp, where J is chosen 
as one of the objectives in eqs. (4)-(6). Evolutionary 
algorithms form the base of many multi-objective 
optimization codes, and here a multi-objective vari¬ 
ant of the covariance matrix adaptation method 
(MO-CMA) was implemented with a population size 
of 100 members attempting to locate the invariant 
front {i.e., the highest flat value of Kp) for each 
type of control objective [8, 26]. Such calculations 
typically required 1000 generations for two-level ob¬ 
jectives, and 10 4 — 10 5 generations for four-level ob¬ 
jectives. The MO-CMA controls were amplitudes a*, 
and phases (f>k of a Fourier field structure: 

£i{t) = '^2a k sin{ujkt + (f>k)- (9) 

fc 

The frequencies were set at uJk £ [0,107r] (2-level 
system) and tok £ [0,207r] (4-level system), and 
A ujk = 7r. The initial amplitudes were randomly 
chosen from either the low fluence regime {ak £ 
[0,0.05]) or a potentially high fluence regime (afc £ 
[0,50]), with fluence defined as Kf i = e 2 (t)dt. 
As an alternative to the multi-objective search al¬ 
gorithm, MC solutions were generated from an en¬ 
semble of D-MORPH optimization trajectories, per¬ 
formed using random initial fields defined by eq. (9). 
In the D-MORPH optimizations, field values e{t) at 
the designated time points were used as controls, in 
contrast to MO-CMA where amplitudes and phases 


3 


are used as controls. Parametrized controls are a 
practical constraint when using MO-CMA, and may 
impact the maximum Kp value. Each MC run re¬ 
sults in a flat, nearly invariant value for Kp (Fig. 
1), and the highest such line is taken as the optimal 
value of Kp drawn from a large set of MC runs. 



FIG. 2: Color online. MO-CMA and Monte-Carlo 
(MC) fronts for dual robustness and fidelity 
objectives, a), c) field noise and b), d) detuning 
noise of two- and four-level objectives, respectively. 
Threshold points Ej = —Kp shown as dots for 
converged results, beyond which noise effects 
dictate that the fidelity J cannot be enhanced 
without first re-engineering the system to improve 
its noise characteristics. Dotted line in c) shows 
MO-CMA performance for comparison with the 
MC outcome. 

Figure 2 shows the best outcomes between MO- 
CMA and MC fronts for optimization of Ej and Kp. 
While a degree of competition is initially present 
at higher Ej, convergence toward invariance is dis¬ 
played at lower Ej between the two costs for any 
objective. In common multi-objective optimization 
problems, MC is typically not an effective algorithm, 
and it is quite surprising that MC solutions out¬ 
performed MO-CMA in several cases, needing only 
modest numbers of MC solutions (Fig. 1 inset). This 
behavior indicates that the top of the quantum con¬ 
trol landscape is often replete with robust solutions 
readily located by moderate number of MC runs. In 
a case where MC performed quite well (Fig. 2 c), 
the MO-CMA front is also shown as a dotted line 
for comparison. For MC solutions in the invariant 
regions, low fluence solutions were common for field 
noise, and high fluence solutions common for detun¬ 


ing noise (not shown). The MO-CMA code located 
fields of high fluence in both cases. 

As mentioned earlier, the implications for the ap¬ 
pearance of an invariance plateau suggest that be¬ 
yond a threshold point Ej = —Kp depicted by the 
dots in Fig. 2, there is no benefit in attempting to 
locate controls that parametrize higher fidelity val¬ 
ues of J, as the effects of noise would overcome those 
efforts. The location of this threshold depends upon 
the objective being optimized, as well as the strength 
and nature of the noise. Attempts to go beyond 
this threshold require re-engineering of the system 
to have a more favorable noise correlation function 
such that the threshold can be lowered to allow for 
a higher fidelity primary objective value for J, as 
demonstrated in Fig. 1. In a more general context, 
any primary objective can only be improved beyond 
its plateaued, invariant secondary threshold value 
through some suitable form of system engineering. 

The invariance behavior between the primary and 
secondary objective still leaves open likely complex, 
competitive behavior among the secondary objec¬ 
tives. A MC sampling of such behavior (ie., the 
bold front) is shown in Figure 3 for Ep H: K e , and 
fluence Kj costs. The invariant behavior will exist 
for any number of secondary objectives, illustrated 
here for Ep H vs. K e or / (shown as projected black 
lines) at all values of the secondary objectives at high 
fidelity of the primary objective Ep H . Interestingly, 
fields that are optimal for Fh and K t are not opti¬ 
mal for Fp and /, indicated by the fold of the 3D 
front at low /. The general features of the 3D front 
are only approximately identified in Fig. 3 through 
modest MC sampling, and including knowledge of 
the landscape principle within the sophistication of 
evolutionary codes could possibly improve search ef¬ 
ficiency for many secondary objectives. 

To summarize, when multi-objective optimiza¬ 
tions focus in the regime of high fidelity for a partic¬ 
ular primary objective, then the landscape principle 
assures flat behavior for all secondary objectives that 
only depend on the control field. In regards to op¬ 
timizing robustness alongside fidelity, observing in¬ 
variance behavior reveals the attractive absence of a 
major obstacle for advancing quantum technologies, 
when it is sufficient to operate at Ej > —Kp. If 
one needs to go beyond the latter point in fidelity, 
an emphasis must be placed on appropriate further 
system engineering. The same conclusion would ap¬ 
ply more generally to encroachment of any invari¬ 
ant secondary objective where the fidelity passes a 
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FIG. 3: Color online. 3D Pareto surface identified 
using MC trajectories (solid blue lines) for Ep H ,K e , 
and / costs. Each MC trajectory shows invariance 
of a second objective (i.e., projected black lines), 
but displays complex 2D behavior between the two 
ancillary objectives (bold line). 


practical threshold of undesirable competition with a 
secondary objective. This paper dealt with systems 
undergoing unitary dynamics, and the same princi¬ 
ple would apply in cases where at least the signifi¬ 
cant portion of the local environment is also under 
full control. Including K u in this work is an effec¬ 
tive stand-in for that situation. Finally, the invari¬ 
ance phenomenon explained here is not restricted to 
quantum dynamics, as the same landscape princi¬ 
ple is also known to be operative for closed, classical 
Hamiltonian systems [27]. Thus, the principle estab¬ 
lished in this paper should have wider applicability 
in multi-objective optimization where a primary ob¬ 
jective is demanded to have high, reachable fidelity, 
and can be exploited in design simulations, as well 
as directly in the laboratory in an iterative fashion. 
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